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! Magnetic elements on the quiet sun are buffeted by convective ffows that 

cause lateral motions on timescales of minutes. The magnetic elements can be 
^ i observed as bright points (BPs) in the G band at 4305 A. We present observations 

^ ! of BPs based on a long sequence of G-band images recorded with the Dutch 

Open Telescope (DOT) and post-processed using speckle masking techniques. 
From these images we measured the proper motions of isolated BPs and derived 
the auto-correlation function of their velocity relative to the solar granulation 
pattern. The accuracy of BP position measurements is estimated to be less than 
23 km on the Sun. The rms velocity of BPs (corrected for measurement errors) 
is about 0.89 km s~^, and the correlation time of BP motions is about 60 s. This 
rms velocity is about 3 times the velocity measured using cork tracking, almost 
certainly due to the fact that isolated BPs move more rapidly than clusters of 
BPs. We also searched for evidence of vorticity in the motions of G band BPs. 



Subject headings: Sun: atmospheric motions — Sun: magnetic fields — Sun: 
faculae, plages — techniques: high angular resolution 
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1. Introduction 

Observations of the Sun with high spatial resolution show network bright points (Muller 
1983, 1985, 1994; Muller & Keil 1983; Muller & Roudier 1984, 1992) and "filigree" (Dunn & 
Zirker 1973; Mehltretter 1974; Berger et al 1995), which are chains of bright features located 
within the intergranular lanes. The bright points and filigree are seen in the wings of strong 
spectral lines such as Ha and Ca II H & K, in lines formed in the photosphere, and even 
at continuum wavelengths (with reduced contrast). The widths of these structures is 100 
to 200 km, at the limit of resolution of ground-based solar telescopes. In the following we 
collectively refer to these bright structures as bright points (BPs). The BPs are associated 
with regions of strong magnetic field (Chapman &; Sheeley 1968; Title, Tarbell & Topka 
1987; Simon et al 1988; Title et al 1989, 1992; Keller 1992) and correspond to magnetic fiux 
tubes of kilogauss field strength that stand nearly vertically in the solar atmosphere (Stenfio 
1973; Stenflo & Harvey 1985; Sanchez Almeida & Martinez Pillet 1994; see review by Solanki 
1993). 

The dynamical behavior of BPs has been studied by a number of authors. Muller 
(1983) found that facular points on the quiet sun are predominantly located in patches at 
the periphery of supergranule cells, indicating that the magnetic elements are advected by 
the supergranular flow. The BPs always flrst appear in the dark spaces at the junction of 
several granules, never inside a granule nor in the space between only two granules. As the 
granulation pattern evolves, the BPs remain in the intergranular spaces throughout their 
lifetime, but not necessarily at the junction of several granules like at the time of their first 
appearance. New BPs have a tendency to appear adjacent to existing points, and 15 % of 
the BPs seem to split into two points which move apart until a separation of 1 to 1.5" is 
reached. 

Muller et al (1994) measured velocities of 29 isolated BPs and found a mean speed of 1.33 
km s^^ Strous (1994) studied BPs in a growing active region. Using line-center images taken 
in Fe I (5576 A), he found velocities between 0.26 km s~^ and 0.62 km s^^. Berger & Title 
(1996) measured velocities of 1 to 5 km s~^ for G-band (4305 A) BPs in the "moat" around 
a sunspot; they showed that the motions are constrained to the intergranular lanes and 
are primarily driven by the evolution of the granulation pattern. They found that the BPs 
continually split up and merge, with a mean time between splitting events of few hundred 
seconds. Lofdahl et al (1998) and Berger et al (1998) analyzed G-band and continuum 
images obtained at the Swedish Vacuum Solar Telescope (SVST) on La Palma and found 
rapid splitting and merging of BPs in an enhanced network region. Van Ballegooijen et al 
(1998) used a "cork tracking" method to measure the proper motions of G-band BPs and 
found that BPs appear to be passively advected by the granulation fiow. The filigree are 
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known to be associated with abnormal granulation patterns (Dunn & Zirker 1973), and the 
granules near network BPs are smaller and more numerous than near a normal intergranular 
space (Muller, Hulot & Roudier 1989), suggesting that the magnetic field has some efi^ect on 
the granulation flow. 

Little is known about the small-scale dynamics of flows in intergranular lanes and the 

interaction of magnetic elements with such flows. Theoretical models (e.g., Stein & Nordlund 
2000; Emonct & Cattaneo 2001) predict that vorticity is concentrated in the lanes and 
that magnetic elements exhibit rapid rotational motions. At present, there is no direct 
observational evidence for vorticity in intergranular lanes. More detailed measurements of 
vorticity are needed for testing magneto-convection models, and as input for models of flux- 
tube waves that heat the upper solar atmosphere (Hasan et al 2002). 

In this paper we present measurements of BP proper motions within the intergranular 
lanes. In particular, we look for evidence of vortical motions. Tracking BPs is a difficult 
task. Long sequences of images with very good seeing (or corrections of seeing) are needed. 
Image jitter due to residual atmospheric effects needs to be reduced as much as possible. The 
measurements of BP positions need to be made with respect to some frame of reference. The 
solar granulation is used as a reference, and the motion of this reference frame is determined 
by correlation tracking on the granulation pattern. We also tested measuring the BPs with 
respect to the average position of the bright points in each field, but found that the result 
is nearly identical to the result from correlation tracking. BPs vary widely in lifetime and 
contrast, so it is very difficult to automate the process of tracking them. In this paper, we 
use manual selection of BPs followed by fitting the peaks. 

2. Observations 

Image sequences recorded with good seeing in the G band (4305 A) are essential for 
tracking BPs over many frames. Here we use G-band images recorded with the Dutch Open 
Telescope (DOT) (Hammerschlag & Bettonvil 1998; Rutten et al 2001, 2002; Siitterlin et al 
2001). The image sequence was collected on 2001 October 19 starting at 10:39 UT under 
good seeing conditions from a Network region close to disk center. The pixel size for these 
images was 0.071" (corresponding to 51.5 km on the Sun), and the filter bandpass was 10 A. 
The data were recorded at a rate of 6 frames per second, and a burst of 100 frames was 
collected every 30 s. Each burst is combined into a single reconstructed frame using speckle 
image processing. The speckle processing uses an improved version of the code of de Boer, 
Kncer & Nesis (1992), who first applied the speckle masking method of Weigelt (1977) to 
solar images, making use of the spectral line ratio technique (von der Liihe 1984). The 
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resulting 199 frame sequence (1 frame every 30 seconds) achieves 0.2" resolution over the 
entire sequence. This is very close to the diffraction limit of the 43.8 cm aperture. Figure 
?? is a typical frame from the sequence. 

The images in the sequence were crudely co-aligned by visual tracking of a group of BPs. 
We then selected two subregions, each 256x256 pixels in size, for measurement of BP posi- 
tions (see Figure ??). For each siibrcgion we used correlation tracking to determine a series 
of displacements of the granulation pattern in consecutive frames. These displacements were 
measured with sub-pixel accuracy by finding the peak of the cross-correlation of consecutive 
subimagcs and fitting a quadratic function to the values near the peak. We accumulated 
these displacements to obtain a series of positions of the granulation pattern, (xgran.n, 2/gran,n), 
where n is the frame index. We then went through each frame many times, using a cursor 
to select a BP and follow it for as many frames as we could, before it disappeared or merged 
with other BPs. Each BP position was again measured with sub-pixel accuracy by quadratic 
fitting, and the granulation position was subtracted. The result was a list of BP positions 
Xk,n and yk,n in a reference frame tied to the solar granulation pattern, where k is the BP 
index. The advantage of this method is that an accurate co-alignment and interpolation of 
the subimages is not required. 

We only selected BPs that were clearly separated from other BPs so that we could follow 
them without confusing their identity. We stopped tracking them if they started to merge 
with other BPs. Any BPs that we could not track for at least 3.5 minute (7 frames) were 
discarded. Our measurements yielded 2992 positions for 161 BPs. The longest-living BP we 
were able to track lived for 25 minutes (50 frames). The mean lifetime was 9.2 minutes. 

Figure 2 shows four subregion frames, separated in time by 1 minute. The circles show 
the positions of the BPs selected from each frame. Figure 3 shows the paths of four individual 
BPs. In some cases, like the upper-left and upper-right plots, the BP followed a simple curved 
path. In other cases, the motions were erratic and complicated, as in the lower- right plot. 
Some BPs move in nearly straight lines, following intergranular lanes. 

3. Analysis 

We calculated BP velocities Vx^k,n and Vy^k,n from position differences over time intervals 
of 30 seconds. Figure 4 shows histograms of v^, Vy and v (= -^v"^ + Vy), and a scatter plot of 
Vx versus Vy. The histograms are consistent with a Gaussian distribution with rms velocity 
of 1.31 km s~^. This includes both real signal and measurement errors. We also computed 
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the auto-correlations Cxx,m and Cyy^m-, and the cross-correlation Cxy,m' 

Cxx,ra i,'^x,k,nVx,k,n+m) 1 Cyy^jn (Vy^k,n'^y,k,n+m) i Cxy,m {l^x,k,n'^y,k,n+m) j (1) 

where m is the number of time steps separating two velocity measurements for the same BP, 
and the brackets denote an average over the BP index k and frame number n. Figures 5a 
and 5b show the auto-correlations as function of the delay time t = 30m [s]. Figure 5c shows 
the cross-correlation function. The velocities Vx and Vy appear to be uncorrclatcd to within 
the measurement errors (note the expanded vertical scale of this figure compared to Figures 
5a and 5b). Figure 5d shows the number of measurements used to derive these correlation 
functions. 

Each auto-correlation function shows a strong peak at zero time lag, and dips at lags of 
±1 frame (30 seconds time delay). We suggest that these features are due to the measurement 
errors. Let Ux,k,n be the x-component of the true velocity of BP k in the time interval between 
frames n and n + 1, then the measured velocity is 

'^x,k,n — "^XjkjU ~i ~/\ f ' ^ ' 

where At = 30 [s] is the time between frames, and 5xk,n are the errors in position mea- 
surement. Now assume that these errors are randomly distributed with standard deviation 
a (the same for all BPs) and are uncorreleated from frame to frame. Then the observed 
auto-correlation Cxx,m is: 

Cxxfi = (ul^k^n) + '^v' 
Cxx,±l = {Ux,k,nUx,k,n±l) — ^v/^j (3) 
Cxx,m {'^x,k,n'^x,k,n+m} fo^ I'^l ^ 2, 

and similar for Cyy^m^ where = a\/2/ At is the error in the velocity measurement. The 
position errors increase the observed correlation at m = by cr^ and reduce the correlation 
at m = ±1 by —(7^/2. Therefore, the average auto-correlations, defined by 

1+1 ~ 1 +^ 

^xx = o ^ V Cxx,m-i Cyy = — ^ ] Cyy^rn-i (4) 



arc not systematically affected by the measurement errors. These average correlations arc 
indicated by the thin lines in Figures 5a and 5b. The difference between Cxxfi and Cxx 
provides an upper limit on cr^ because the true auto-correlation {ux^k,nUx,k,n+m) decreases 
with increasing |m|. The auto-correlations for Vx and Vy yield upper limits on (jy of 0.98 km 
s~^ and 1.12 km s~^, respectively. Therefore, the upper limit on the position error a is about 
24 km. 



-6- 



Figures 5a and 5b show that the auto-correlations are significantly different from zero 
for delays up to 500 s. The average correlations for are C^x = 0.83 ± 0.03 km^s~^ and 
Cyy = 0.78 ± 0.03 km^s~^, where the error estimates are based on the cross-correlations 
shown in Figure 5c. This provides a lower limit on the true velocity dispersion Uq of BPs in 
one direction: Uo ~ 0.89 km s~^ based on the average of Cxx and Cyy. Furthermore, Cxx and 
Cyy are about twice the auto-correlations at m = ±2. Therefore, the correlation time of the 
BP velocity is about 60 seconds (2 frames). 

We were interested in using BPs as a possible diagnostic of vortical motions in the 
intergranular lanes. A BP caught up in a vortex within an intergranular lane might exhibit 
rotational motions. In an effort to detect such motions, we computed changes in direction 
angle, A0, of BP velocity in two neighboring time intervals. For this calculation we discarded 
any measurements for which v < 1.5 km s~^ in either time interval, since the angle errors 
became large for smaller velocities. This left us with 687 measurements. A histogram of 
is shown in Figure 6a. The distribution shows a broad peak at A0 = 0, indicating that the 
direction of motion persists from one time interval to the next. Figure 6b shows a histogram 
of the centrifugal acceleration, vd(j)/dt, which is the relevant quantity when considering the 
generation of transverse waves in magnetic flux tubes. Histograms such as these could be 
useful tools for comparing models of flux tube waves to observations. 

The above measurements use the granulation pattern as a reference frame. In order 
to check whether the granulation provides a stable reference frame, we also measured BP 
positions relative to each other (independent of the granulation pattern). For this purpose 
we computed a reference position {xrei,n, Uvci-n) as follows. For each pair of frames n and 
n -|- 1, we determine the average positions {xn,yn) and {xn+i,yn+i) of all BPs in common 
between these consecutive frames. Then we accumulate the position differences: 

•^ref,n+l — -^refjn ~l~ -^n+l -^m Z/ref,n+l — yTei,n ~l~ Vn+l Vm (5) 

starting with Xref,o = Z/ref,o = in the first frame. This provides a reference position against 
which the individual BPs can be measured. The relative positions are a;'^,^ = Xk,n ~ x^ef,n 
and y'l^n = yk,n — J/rcf.n- The auto-correlations of v'^ and v'y yield C^x = 0.71 ± 0.03 km^s~^ 
and Cyy = 0.66 ± 0.03 km^s^^, which is about 15% less than the values derived for the 
granulation reference frame. This reduction in BP velocities is to be expected because there 
are on average only 7.0 BPs in common between consecutive frames, and subtracting the 
mean velocity of these 7 BPs should reduce Cxx and Cyy by a factor 6/7. The fact that 
the observed reduction is not any larger than predicted implies that there is nothing to be 
gained by measuring the BP motions relative to each other, and in fact this causes the BP 
velocities to be systematically underestimated. We conclude that the granulation pattern 
provides a stable reference frame for measurement of BP motions. 
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4. Discussion 

We have shown that with good enough data we can precisely measure positions for 
isolated BPs, obtaining both their linear and angular motions. The errors due to image 
jitter have been reduced to levels of 25 km (0.5 pixels) or less. This allows us to measure 
the motions of individual BPs to a precision of less than 1 km s~^ for each 30 second time 
interval, which will be important for studying the generation of MHD waves in flux tubes. 
We also track BPs for sufficiently long times that we can see the regularity in their motions. 

The rms velocities we measured for isolated BPs were almost a factor of three higher 
than velocities measured using cork tracking (van Ballegooijen et al 1998). We beheve this 
is due to a selection effect. By only measuring isolated BPs we measure the fastest-moving 
features. Once BPs cluster with other BPs, mostly at lane interstices, they seem to be 
"captured" by the group and their motions are reduced. There are many more BPs in these 
clusters than there are isolated BPs so the cluster statistics dominate when measurements 
are made with cork tracking and the average velocity is reduced. 

In measuring the motions of BPs, we have found that they move in somewhat circular 
paths. Combining the angular change of their motions with the distance they travel is a 
potential way of estimating the (vertical) vorticity in the flow field, assuming the BPs act 
as test particles. While performing the tracking of individual BPs, we saw little evidence 
of vorticity other than the tracks of the BPs themselves. We looked for, but did not find, 
pairs of BPs that orbited one another; in fact, BPs that came in close proximity to other 
BPs tended to have reduced motion. We also saw little evidence of swirling motions in the 
granules or granular lanes that would correspond to the observed BP motions. This might 
indicate that BPs are affected not only by surface flows associated with the solar granulation, 
but also by other flows that occur at larger depth below the photosphere. Further modeling 
will be required to establish whether or not the BPs move passively with the surface flow, 
and to determine whether the observed BP motions are consistent with models of flux tube 
waves and heating of the upper solar atmosphere. 

The authors wish to thank R. J. Rutten, R. Hammerschlag and F. Bettonvil for the use 
of DOT solar data. We also wish to especially thank R.J. Rutten for his very helpful sugges- 
tions on an earlier draft of this paper. We wish to acknowledge support for the work from 
the National Science Foundation (Grant ATM-9811523) to the Smithsonian Astrophysical 
Observatory. 
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Fig. 1. — G-band frame taken at 10:39:35 UT. The boxes show the two subregions used for 
measurement of BP positions. 
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Fig. 2. — Four images and with BP locations circled. The frames are separated in time by 
1 minute. 
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Fig. 3. — Four bright point tracks. The starting point of each track is indicated by a star. 
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Fig. 4. — Histograms of measured BP velocities: (a) v^; (b) w^; (c) v. The results are 
compared with a Gaussian distribution with rms velocity 1.31 km s~^. Panel (d) shows 
plotted against Vy. 
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Fig. 5. — Correlation functions of BP velocity and Vy. (a) Observed auto-correlation Cxx,m 
as function of delay time t = 30m [s] {thick curve), and the average auto-correlation C 
{thin curve), (b) Similar for the observed auto-correlation Cyy^m {thick curve) and average 
correlation Cyy {thin curve), (c) Cross-correlation C^y^m as function of delay time t. (d) 
Number of measurements per bin used in panels (a), (b) and (c). 
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Fig. 6. — (a) Histogram of the change A0 in direction angle of BP velocity, (b) Histogram 
of the centrifugal acceleration vd(j)/dt. 



